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The energy dependence of hadron production in central Pb+Pb collisions is presented 
and discussed. In particular, midrapidity m^-spectra for n~ , K~, K + , p, p, d, 4>, A and A 
at 40, 80 and 158 AGeV are shown. In addition H and Q spectra are available at 158 AGeV. 
The spectra allow to determine the thermal freeze-out temperature T and the transverse 
flow velocity (3t at the three energies. We do not observe a significant energy dependence 
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of these parameters; furthermore there is no indication of early thermal freeze-out of H 
and Q at 158 AGeV. Rapidity spectra for ir~, K~, K + and (f) at 40, 80 and 158 AGeV 
are shown, as well as first results on Q rapidity distributions at 158 AGeV. The chemical 
freeze-out parameters T and at the three energies are determined from the total 
yields. The parameters are close to the expected phase boundary in the SPS energy range 
and above. Using the total yields of kaons and lambdas, the energy dependence of the 
strangeness to pion ratio is discussed. A maximum in this ratio is found at 40 AGeV. This 
maximum could indicate the formation of deconfined matter at energies above 40 AGeV. 
A search for open charm in a large sample of 158 AGeV events is presented. No signal is 
observed. This result is compared to several model predictions. 

1. Introduction 

In the last few years NA49 has taken data at 40, 80 and 158 AGeV beam energy to study 
the energy dependence of particle production in nucleus-nucleus collisions and to search 
for indications of the onset of deconfinement in these collisions. In this paper recent results 
on particle spectra and yields are reported. Results by NA49 on other observables such 
as HBT, fluctuations and system size dependences from p+p, p+A, and A+A (peripheral 
Pb+Pb, central Si+Si and C+C) are presented elsewhere in these proceedings jll'2l,'j|4j . 

In addition to the energy scan programme, a large dataset of central Pb+Pb data has 
been collected at 158 AGeV, allowing to study rare particles such as Q, D, and </> — > e + e~. 
Results on Q and open charm production are presented here. 

2. Experiment and centrality selection 

The NA49 experimental setup jH] consists of four TPCs, two of which are in a magnetic 
field. The Main TPCs (MTPC) are outside the magnetic field and are used to identify 
particles by measuring the energy loss dE/dx. The particle identification capabilities of 
the MTPCs are augmented by two time-of-flight (TOF) detector arrays which have accep- 
tance around midrapidity for kaons and close to midrapidity for protons and deuterons. 
Centrality selection of the events is done using a calorimeter which detects the projectile 
spectators. At 40 and 80 AGeV the 7.2% most central events were selected. For the 
158 AGeV data the on-line centrality cut is at 10%, but for the analysis of pions, kaons 
and (p an off-line cut was made at 5%. The high statistics data sample used for the Q 
analysis at 158 AGeV was taken with a centrality trigger at 20% of the inelastic cross 
section. 

3. Transverse mass spectra and thermal freeze-out 

In Fig. ^ all available results on m^-spectra as measured in central Pb+Pb collisions 
by NA49 are collected. 2 The spectra are measured at or close to midrapidity except for 
the 0, H and Q which are integrated over approximately 1 unit of rapidity, to increase 
statistics. 

2 The p, Sl~ and £l + mr-spectra at 158 AGeV appearing in the printed version of these proceedings (Nucl. 
Phys. A715 161c) are unfortunately incorrect and therefore differ slightly from those presented here. 
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Figure 1. Transverse mass spectra in central 40 (left), 80 (middle) and 158 AGeV (right) 
Pb+Pb collisions. The lines show the result of the transverse flow fit as described in the 
text. 
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The H and <fi at 158 AGeV have already been published |6|7j while the kaon and pion 
spectra are submitted for publication [H] . All other presented results are still preliminary. 

In a hydro-dynamical picture, the mr-spectra of particles are sensitive to transverse 
flow. To characterise the flow, the spectra were fitted with [H] 

dN frriT cosh p\ f p T smhp\ 

oc m T K x J . (1) 



mTdrriTdy \ T J y T 

A combined fit of several spectra with this function allows to determine the thermal freeze- 
out temperature T and the mean transverse flow velocity (3t (p = atanh ftp)- The ti~ 
and deuteron spectra were excluded from the fit (dashed lines in Fig. because the pions 
are expected to have a significant contribution from resonance decays and the deuterons 
may be formed by coalescence. The particles (baryons and K + ) and anti-particles (anti- 
baryons and K~) were fitted separately. The was included in the anti-particle fit, 
because it is not expected to be sensitive to the baryon density. The fit describes the data 
reasonably well, even including the pions at higher my and the deuterons. The relatively 
large \ 2 P er degree of freedom can partly be explained by the fact that only statistical 
errors were taken into account in the fit (systematic errors are not yet available for all 
spectra). Since the slopes of the m T -spectra of particles and anti-particles are similar, the 
differences between the freeze-out parameters for both groups of particles are small. 

The differences between the three different energies are not larger than the difference 
between the particles and anti-particles at each energy, implying that the energy depen- 
dence of the thermal freeze-out conditions is small at these energies. 

The 5 and Q spectra, which are only available at 158 AGeV, do not show large devia- 
tions from the fit. Thus, based on these data, there is no indication of an early freeze-out 
as has been suggested by several authors |10|11| . The slope parameters of the Q spectra 
are compatible with earlier measurements by WA97 [T2"] . 



4. Rapidity spectra and baryon density 

The m^-spectra of kaons as shown in Fig. ^ were obtained using combined dE/dx and 
TOF information at midrapidity. The dE/dx measurement in the TPC alone allows to 
identify kaons at forward rapidities. The resulting spectra are shown in Fig. together 
with 7r~ distributions which were obtained using unidentified negatively charged particles, 
corrected for contributions of K~,p and non- vertex tracks. 

The K~ distributions are narrower than the K + distributions at all energies. This could 
indicate that the kaons are sensitive to the local baryon density, which is reflected in the 
difference between the A and A spectra as well (for details on the A analysis, see ^Hj in 
these proceedings). 

The spectra are obtained from invariant mass distributions of K + K~ pairs. The 
resulting rapidity spectra are shown in Fig. El where the new results at 40 and 80 AGeV are 
compared to the published result at 158 AGeV |7| . Clearly, the width of the distribution, 
as well as the total yield, increases with beam energy. 

The large statistics data sample of the 20% most central Pb+Pb collisions at 158 AGeV 
beam energy, which was taken in the year 2000, has been used for the Q analysis. The 
resulting m^-spectra are shown in Fig. ^ and the corresponding rapidity spectra in Fig.0] 
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Figure 2. Rapidity spectra for kaons and pions at three different energies. The open 
points are reflected around midrapidity. For kaons, the circles indicate the result of the 
dE/dx analysis, while the squares are the result of the combined dE / dx-TOY analysis. 
The distributions are fitted with a double Gaussian and the resulting full width at half 
maximum (FWHM) is indicated in the figure. 
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Figure 3. Rapidity spectra of production at 40, 80 and 158 AGeV respectively. The 
open points are reflected around mid-rapidity. The full lines show a Gaussian fit, used to 
obtain the mean multiplicities in full phase space (<fi). The fitted widths a are given as 
well. 
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Figure 4. Rapidity spectra of fl~ and Q + at 
158 AGeV. The total yields (fl) were calculated by in- 
tegrating the fitted Gaussians (full line). The widths 
a of the fitted curves are also given. 



Figure 5. Chemical freeze-out 
points as determined from a 
hadron gas fit at different ener- 
gies [To] . The crossover curve is 
taken from |16j . 



Integrating over rapidity gives mean multiplicities per event of = 0.47 ± 0.07 and 
= 0.15 ± 0.02. The observation that the ratio Cl + /Or is below unity cannot be 
explained in a pure string fragmentation model but fits to expectations in an equilibrium 
hadron gas model ^3]. Also note that the Q~ has a slightly wider rapidity distribution 
than the This indicates that even the Q is sensitive to the effect of baryon density. 

5. Total yields and chemical freeze-out 

The measured rapidity spectra have been integrated to obtain the mean multiplicities in 
full phase space. To determine the chemical freeze-out parameters at the three different 
energies, the total yields were fitted with a hadron gas model with partial strangeness 
saturation The obtained temperatures for chemical freeze-out (145-160 MeV) are 
larger than the thermal freeze-out temperatures (115-140 MeV). The chemical freeze-out 
temperature T and baryon chemical potential \xb are compared to results at AGS and 
RHIC in Fig. 03 The freeze-out points follow a smooth curve of increasing temperature 
and decreasing baryon density with energy. The observed freeze-out parameters are close 
to the calculated phase boundary {16J at SPS energies and above. 

6. Energy dependence of strangeness production 

The mean multiplicities of the most abundant carriers of strangeness, kaons and lamb- 
das, are plotted as a function of the centre of mass energy in Fig. |H1 Only the charged 
kaons are shown because the neutral kaons are not measured at all energies. All yields 
are divided by the total pion yield (n) = 1.5((7r + ) + It is seen that the (K + ) / (n) 

ratio peaks at 40 AGeV, whereas (K~) / (tt) increases monotonically. On the other hand, 
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Figure 6. Total yields of K + and A (left) and K~ and A (right), normalised to the 
total pion multiplicity, as a function of y/s for central Pb+Pb (NA49) and Au+Au (AGS) 
collisions. 
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Figure 7. Energy dependence of (K + ) / (n) 
compared to different model predictions 
which do not invoke QGP formation (see 
text). 



Figure 8. Comparison of the strangeness to 
entropy ratio Eg to a model calculation [20J 
which does assume that a QGP is formed 
above approximately 30 AGeV. 
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the A production also peaks around 40 AGeV, whereas the A yield increases monotoni- 
cally with energy. Apparently, both the total s production and the s production have a 
maximum close to 40 AGeV. The maximum in the s production is directly seen K + /tc 
ratio because the contribution of A is small. Although the K~ represent a large fraction 
of the s production, the maximum in the s production is reflected in the A/tt ratio. 

In Fig. [7|the energy dependence of the (K + ) / (tt + ) ratio is compared to several models 
which do not assume the formation of a QGP in the collision. The UrQMD model ^\ 
fails to describe the data because the pion yields are over-estimated by about 30%. The 
RQMD model [TS] gives a much better descriptions, but over-estimates the ratios at 80 
and 158 AGeV. The same is true for the hadron gas model [T^j. In this model, the 
(K + ) / ratio peaks around 40 AGeV due to the interplay of the decreasing baryon 
density and the increasing temperature with energy. 

In Fig. El the energy dependence (F = (-v/i — 2mAr) 3 / 4 / v /i 1//4 ) of the strangeness to 
entropy ratio E s = ((A) + (K) + (K)) / (n) is compared to a model |2Dj which assumes 
that above a certain energy a QGP is formed in the collision. The peak in the model 
curve is caused by this transition. In [5| it is shown that the model also describes the 
energy dependence of pion production. 

NA49 intends to complete the energy scan by taking data at 20 and 30 AGeV in 2002. 

7. Open charm upper limit 

All available central Pb+Pb data at 158 AGeV (900k events at 10% and 3M events 
at 20% centrality) were used in an invariant mass analysis to measure open charm. The 




m(K ,7t + ) (GeV) m(K, jt) (GeV) 

Figure 9. Uncorrected invariant mass distribution of all K~,7r + pairs (left) and invariant 
mass distribution for the sum of all D° and D° candidates after background subtraction 
(right), corrected for acceptance and efficiency The curves indicate predictions of various 
models (see text). 
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invariant mass was calculated for pairs of pion and kaon candidates, to detect the decays 
D° -> K~ + 7T+ and D° -> + tT, 

which have a branching fraction of 3.83%. For tracks which reach the MTPC, the kaon 
candidates were selected by dE/dx measured in the MTPC. For tracks which do not 
reach the MTPC, all candidates were used. The invariant mass histogram for D° before 
background subtraction is shown in the left panel of Fig. El The result after adding 
the D°, subtracting the background and correcting for acceptance and efficiency losses is 
shown in the right panel of Fig. El No signal of the decays is visible. Detector simulations 
have shown that the shape of the expected signal can be parametrised as a Breit-Wigner 
with a width T = 11 MeV. Using this parametrisation, predictions by different models 
are shown in the figure. It is seen that thermal production of charm quarks in a Quark 
Gluon Plasma at T = 265 MeV, as described in [2T], is excluded (dotted curve in Fig. EJ). 
The expectation from the quark coalescence model ALCOR [221 is at the limit of the 
sensitivity of the present analysis (dashed curve). Lower values of the charm yield, such 
as the expectation from perturbative QCD, or a factor 3 enhancement j2H| (full curve) 
cannot be excluded. In order to quantify the result, the data have been fitted with the 
parametrised peak shape, leaving the normalisation as the only free parameter. The result 
is (D° + D°) = -1.07 ±0.84. 

8. Summary and outlook 

Spectra of n, K,p, d, <p, and A in central Pb+Pb collisions at 40, 80 and 158 AGeV are 
presented, together with new results on Q rapidity and vht spectra at 158 AGeV. 

The m^-spectra are compatible with radial flow, with similar temperature T m 130 MeV 
and transverse flow velocity (3t ~ 0.45 at the different energies. The S and Q spectra, 
which are only measured at 158 AGeV, do not show indications of early freeze-out. 

The chemical freeze-out parameters as determined from the measured total yields ap- 
proach the phase boundary in the SPS energy range. 

The energy dependence of strangeness production shows a maximum in the relative 
strangeness yield at 40 AGeV. This maximum is observed in the K + /ir ratio for the s 
and in the A/tt ratio for the s quarks. The maximum can be explained by the onset 
of deconfinement in the energy range between AGS and SPS. Upcoming runs at 20 and 
30 AGeV will provide a more accurate localisation of the maximum. 

To investigate open charm production in heavy ion collisions, an invariant mass analysis 
has been performed using a large sample of 158 AGeV Pb+Pb events. No signal has been 
observed. The result excludes equilibrium charm production in a QGP of T = 265 MeV. 
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